Continuous Crystallization in Hexagonally-Ordered Materials 
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We demonstrate that the phase transition from columnar-hexagonal liquid crystal to hexagonal- 
crystalline solid falls into an unusual universality class, which in three dimensions allows for both 
discontinuous transitions as well as continuous transitions, characterized by a single set of exponents. 
We show by a renormalization group calculation (to first order in e = 4— d) that the critical exponents 
of the continuous transition are precisely those of the XY model, giving rise to a continuous evolution 
of elastic moduli. Although the fixed points of the present model are found to be identical to the 
XY model, the elastic compliance to deformations in the plane of hexagonal order, p, is nonetheless 
shown to critically influence the crystallization transition, with the continuous transition being 
driven to first order by fluctuations as the in-plane response grows weaker, /j-»0. 



On its way from an isotropic fluid to a three- 
dimensional (3D) crystal, a molecular system may pass 
through a nearly limitless sequence of mesophases of 
intermediate symmetry and bulk properties [l|. For 
many systems close to crystallization, a liquid crystalline 
state with only a single continuous symmetry remaining 
appears, the columnar phase, in which molecular con- 
stituents adopt truly long-range two-dimensional (2D) 
order in the plane but remain fluid in the out-of-plane 
direction. This "minimally fluid" state is marginal to 
a full 3D solid in the sense that it is missing only a 
single elastic constant associated with shear along the 
columnar axis. Historically, the columnar-hexagonal 
phase has been studied in the context of discotic, or 
disk-like, molecules [2J, |3|. More recently, an interest 
in the complex ordering of stiff biopolymers, such as 
DNA [4U5 Ja, Lfl] , and filamentous protein aggregates and 
viruses [8|,l9|, has refocused attention onto the specialized 
properties of columnar ordering and its role in biological 
assemblies. Indeed, recent studies predict a direct con- 
nection between the type and quality of hexagonal order- 
ing and the structural [101 ] and mechanical [11,] properties 
of bundle aggregates of bio-filaments, such as those that 
constitute the cytoskeletal networks of living cells. 

In this Letter, we analyze the critical properties of a 
system on the verge of losing the last vestige of fluid 
symmetry. We predict that the columnar-to-solid tran- 
sition may take place either 1) as a second-order transi- 
tion characterized by the critical exponents of the XY 
model or 2) as a fluctuation-driven first-order transi- 
tion, quite analogous to the finite temperature normal- 
to-superconductor transition in metals. The nature of 
the transition is controlled by in-plane elastic response 
of the system, (i, although curiously the coupling to in- 
plane modes of the system is not relevant to linear order 
at the fixed points of the renormalization group flow. Be- 
low we give arguments to suggest that continuous crys- 
tallization can be associated with filamentous structure 
of molecular constituents, while a discontinuous transi- 
tion will be more likely observed for disk- like molecules. 



Finally, we summarize the universal "critical elasticity" 
that emerges as a consequence of a continuous crystal- 
lization transition. 

Dual to the columnar phase is the smectic-A phase of 
rod-like molecules, which has broken translational sym- 
metry along one axis and remains fluid in the 2D plane 
perpendicular to that axis. Owing to the critical im- 
portance of fluctuations of "one-dimensional solids" [l2j , 
the nematic-to-smectic-A transition is notoriously com- 
plex and has received considerable theoretical atten- 
tion [ljj LLJ, UM, [l6| • In comparison, little is known about 
the statistical mechanics columnar-to-solid transition de- 
spite the fact that-as we will demonstrate-the critical 
properties of the former system may be considered a par- 
ticular case of the present, more general, model describ- 
ing the onset of one-dimensional periodic order. Indeed, 
the statistical mechanics of the columnar-to-solid transi- 
tion are surprisingly complex, leading to critical behavior 
with universal and non-universal aspects. 

A columnar liquid crystal is characterized by its elas- 
tic response to mechanical deformations. Because it 
is translationally invariant along the columnar axis, a 
columnar system can be parameterized entirely by a two- 
component displacement field, u(x), which describes the 
in-plane deviation of the column or filament backbones 
from their local equilibrium. Assuming that columns are 
oriented along the z axis, the unit tangent vector repre- 
senting backbone orientation is simply, t = z + d z u. Ow- 
ing to the 2D hexagonal order, columnar systems resist 
in-plane mechanical strains like any isotropic 2D solid. 
Hence, the Hamiltonian describing the elastic response 
of a columnar system has the following form 



ia, 



w± 



d 3 x[A u\, k + 2fi UijUij + Ki(Vj_ ■ t) 2 

+ K 2 (v ± xi) 2 + K 3 (d z i) 2 ]. (1) 



Here, Uy- = (diUj + djUi)/2 is the linearized 2D strain 
tensor, A and fi are the Lame elastic parameters pe- 
nalizing in-plane compression and shear, and K\, K2 
and K3 are the Frank elastic parameters for splay, twist 



and bend deformations, respectively. While in a three- 
dimensional solid one may ignore the higher-order deriva- 
tive contributions represented by the Frank elastic terms, 
in the columnar phase it is necessary to maintain at least 
Kz 7^ as the system is especially soft to bending defor- 
mations (i.e. d±u = but d z u ^ 0). Such a system is 
stable to long-wavelength thermal fluctuations, and both 
orientational order and positional in-plane order are truly 
long-ranged in the thermodynamic limit. 

The crucial distinction between 2D columnar vs. 3D 
crystalline order is the broken translational symmetry 
along the column axis of the latter. We therefore con- 
struct a Landau-Ginzburg free energy in terms of a com- 
plex order parameter, ^>(x), representing broken transla- 
tional order along the columnar direction. Specifically, 
we may define it in terms of the density variation of 
the system near crystallization, 5p(x) ~ Spn cx (x)(l + 
-0(x)e' sz +cc), where 2ir/g is the wavelength of the den- 
sity modulation along the backbone direction. Here, we 
focus on the simplest case, namely the transition from 2D 
hexagonal-columnar order to 3D hexagonal-solid order. 
Clearly, ip(x) represents a smectic-like ordering, with re- 
spective in-plane and out-of-plane gradients of the phase 
of if) corresponding to bending and compressive phonons 
of the nascent ordering. As in smectics, rotational in- 
variance of the system requires a specialized coupling of 
in-plane derivatives of ip to the displacements of column 
orientation. With this in mind, the most general model 
one can write down (to second order in derivatives) to 
describe this transition has the following form, 



1 
H =2 



d 3 x 



rH 2 + -h/f 



+ Cj_|(Vl + igt^ 2 + q, |(0, + tg'u kk )if>\ 2 

+ (SX u\ k + 25fi UijUij) \ip\ 2 . (2) 

With the exception of the terms coupling to the in- 
plane strain tensor, this Hamiltonian has precisely the 
form required to describe the nematic-to-smectic A tran- 
sition [18|. The critical properties of that model have 
been studied extensively, both in the context of liquid 
crystals 1J, [l5| as well as the normal-to-superconducting 
transition in metals [13j, and it was shown that as the 
critical temperature is approached (r —* 0) coupling of 
i\) fluctuations to fluctuations of t drive the transition 
weakly first order. The origin of this effect can be traced 
to the fact for the nematic-to-smectic transitions, fluctu- 
ations of the director field are massless, especially soft, 
as their energy cost goes to zero roughly as fc 2 |tj^(k)| 2 . 
In the present case, only bending modes are massless in 
this sense, as the in-plane ordering has quenched most 
of the long-wavelength fluctuations of the director field. 
Hence, taking /i — > reduces the columnar phase to nc- 
matic so that this limit recovers the first-order nematic- 
to-smectic A transition. On the other hand, in the ab- 



sence of out-of-plane positional order, orientational fluc- 
tuations coupling to order parameter fluctuations go to 
zero as (|t ± | 2 ) - fc s l>- 3 / 4 A 3 ~ 1/4 for large fj, [H]. Intu- 
itively, one would expect that in the limit that in-plane 
fluctuations are "frozen out" by the stiff elastic response 
of the system, the phase transition is restored to second- 
order, as mean-field theory predicts. 

To analyze the critical properties of this system more 
rigorously, we performed the renormalization group (RG) 
analysis of eqs. (Q} and ^ in d = 4 — e dimensions. For 
these purposes it is useful to rescale fields and dimen- 
sions to obtain a more isotropic gauge interaction with 
C|| = C± = C and g' = —g. In this case, fluctuations of 
if> couple to a transverse vector field, n, with n^ = d z u 
and n z = — Vj_ • u. We choose to rescale lengths un- 
der RG transformation by z' — e z and Xj_ = e Xj_, 
where S will be chosen to maintain the isotropic form 
of the gauge coupling. At each RG step, we use stan- 
dard diagrammatic methods, to calculate the perturba- 
tive corrections to the renormalized model from field fluc- 
tuations for large wavevectors Ae~ £ < k < A. The full 
details of this analysis will be presented elsewhere. It is 
nevertheless important to note here that we generically 
focus our analysis on the case when A^A 2 / /i <C 1 (i.e. 
small penetration depth). It is straightforward to show 
from the rescaling procedure above that di ln[A^A 2 / fi] = 
—2 + 4(1 — Si) + 0(e). A careful analysis shows that 
1 — 8i < and generically flows to zero, hence all flows 
tend to the AiA 2 /^ — > limit. 

Although the analysis of this model involves many pa- 
rameters, only a few parameters contribute to the fixed- 
point behavior. Both S/j and SX constitute irrelevant 
parameters under RG flow, and we drop them from the 
subsequent analysis. The most relevant couplings to the 
in-plane elastic response are captured by the parameters, 



V I27r 2 fii2 + a e V h 



(3) 



and Att 2 R^ = gj/ K%. Here, at — X^/fie and we have set 
C( = \. Given these definitions the RG flow is governed 
by the following recursion relations, 
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with dextjdl = 0(e 2 ). Clearly, Qi — > for all flows near 
4 dimensions, and therefore, Re flows to a stable fixed 
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FIG. 1: The RG flows predicted by eq. g]for e = 1, A = a = 

and Rg at its stable fixed point. The blue curves show param- 
eter flow into the stable, Wilson-Fisher fixed point, while the 
red curve shows runaway flow. Dividing stable and unstable 
flow regions is a separatrix (shown in purple) flowing into the 
Gaussian fixed point. The inset show the seperatrix flow for 
both a — and a — » oo far from the fixed point. According to 
these predictions, continuous crystallization is only possible 
for sufficiently small (large) values of Qq (/zo). 



point, R* = 6e. Therefore, the fixed points of eq. ((3]) 
are identical to the isotropic XY model with a Gaussian 
fixed point, r* = u* — 0, and a stable Wilson-Fisher fixed 
point, r* — — eA 2 /5 and u* = 4-7r 2 e/5. Indeed, to linear 
order, the couplings to in-plane elasticity are irrelevant 
at the fixed points, leaving the critical exponents of those 
points unchanged. 

In spite of the irrelevance of the gauge couplings at the 
fixed points, a very small non-linear perturbation, 0(e), 
around the Gaussian fixed point strongly alters the RG 
flow. Focusing on the evolution of the quartic ij> coupling 

3/2 

in the vicinity of R = R*, the coefficient of the Q e term 
is of order e 1 ' 2 , so that it generically dominates the flow 
in a neighborhood of order e of that point. In the limit of 
small (-C e) ug and Qg, the solution to the recursion rela- 
tion for ug shows that RG flows are only stable provided 

3/2 

that uq > eCQ , where uq and Qq are bare parameters 
and C is a non-universal constant, dependent on a. Out- 
side of this region, flows runaway to ug — > — oo, which we 
interpret as a signature of a first-order transition. Hence, 
this analysis predicts regions of stable and unstable RG 
flow for d < 4. In particular, we show numerical solu- 
tions to the flow equations e = 1 in the Qg — ug plane in 
Figure Q] (evaluated for Rg = R*). A separatrix flowing 
into the Gaussian fixed point divides stable and runaway 
flows, indicating that the nature of transition, first or 
second order, is governed by the intial values uq and Qo- 
Interestingly, a similar RG flow is found for zero temper- 
ature, low-dimensional superconducting systems in the 
presence of long-range repulsive interactions [20[. Not 
unlike the effect of in-plane order in columnar systems, 
the effect of long-range interactions is to sufficiently sup- 



press critical gauge fluctuations of the superconducting 
system, allowing a phase transition to proceed as second 
order. 

Interestingly, we see from Fig. [T]that according to this 
analysis in d — 3 stable RG flow is only possible below a 
maximum value of Q, which ranges from 0.08 for a = 
to 0.21 for a — * oo. Hence, this predicts a critical value 
of in-plane shear modulus, /z, that separates first-order 
crystallization from second-order crystallization transi- 
tions. Existing experimental studies of the columnar- 
to-solid transition confirm the existence of both contin- 
uous 21 1 and discontinuous [22| transitions. In order 
to understand what distinguishes the thermodynamics of 
these systems, we estimate the value of in-plane shear of 
a columnar system in terms of the microscopic in-plane, 
a±, and out-of-plane, an, length scales: zx ~ a7 a ± 
(in units of ksT). Here, a\\ ~ 2-7T / ' g is the layering dis- 
tance of smcctic-like order along the column direction and 
a±_ corresponds to the in-plane separation of columns. 
Using the eq. ([3]) and the fact that our analysis is car- 
ried out near d = 4, we estimate that the crystallization 
transition is first-order provided that a\\/a± < (AajJ 2 . 
Based on this estimate, we expect that crystallization of 
high-aspect ratio, filament-like molecules will be continu- 
ous, while crystallization is discontinuous for low-aspect 
ratio, disk-like molecules. Indeed, this picture is in per- 
fect agreement with the respective results of DNA [21[ 
and discotic molecule [22j crystallization studies. 

These results establish that, to order e, the crys- 
tallization of columnar systems may occur either as a 
fluctuation-driven first-order transition or as a continu- 
ous, second-order transition characterized by XY expo- 
nents, i.e. £ ~ |t| _I/ with v ~ 2/3. Here, t is a generalized 
critical parameter, which may reflect either changes in 
temperature [22J or concentration [9|, [2lJ . We now briefly 
discuss the generic consequences of the continuous onset 
of broken translational symmetry in three dimensions, 
focusing on the response of the system to shear deforma- 
tions that tend to slide the columns past one another. 
This issue has been discussed previously in the context 
of a specific crystallization transition occurring in poly- 
electrolyte complexes, in which the elastic response to 
an external strain on the ci-dimensional classical system 
was mapped onto the current response to an external vec- 
tor potential, Aj_, of a zero-temperature superconducting 
system in d — 1 dimensions [23j . 

In the most general case, the distinct properties of 
columnar and solid order are encoded in the elastic re- 
sponse /i z j_(k) to an inhomogeneous strain u^j_(k). Be- 
cause a columnar system is translationally invariant, it 
has no response to uniform (k = 0) strain. In the solid 
phase, the long-wavelengh response to non-uniform dis- 
placements along the z direction has the form, 
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FIG. 2: A sketch of the singular elastic behavior of colum- 
nar systems undergoing continuous crystallization transitions. 
The Frank elastic constants diverge at the critical point, while 
a response to uniform shear grows continuously with t. At 
t — generic scaling arguments predict a non-analytic, criti- 
cal elasticity. 



where <p is the phase of the complex order parameter, 
ijj. By dimensional analysis ^ z j_(k = 0) ~ £-( rf - 2 ) 5 
since critical contributions to the free energy are gov- 
erned entirely by the divergent correlation length. Fol- 
lowing a scaling argument of Fisher and coworkers devel- 
oped for quantum phase transitions of two-dimensional 
superconductors [24j | we may deduce the wavelength de- 
pendence of /^ z j_(k) near the critical point just below 
the crystallization transition. Specifically, we must have 
MzJ-(k) ~ €~ f(k£)i where f(x) is a dimensionless 
function which goes to a constant for x = and must al- 
low /Ltjsj_(k) to remain finite as x — > oo since the response 
of the system to an inhomogenous k ^ shear will be 
finite at the critical point . This requirement is satis- 
fied by f(x) = 1 + Cox^ d ^ 2 \ which predicts an unusual 
critical elastic response to shear /U z j_(k) ~ |k| in three 
dimensions when £ — > oo. Approaching the critical point 
from the columnar phase where <f> has no meaning, the 
singular contribution to the elastic energy of the system 
can be written as, 



sn c 



ilumnar 



d d k 

(2tt)° 



^±(k)|t ± (k)| 



(6) 



Since the long-wavelength response of the columnar 
phase to orientational deformations goes as fc 2 |tj_(k)| 2 , 
we deduce that f(x) oc x 2 above melting and [i z ±(k) ~ 
£-(d-4)^2 ^_ e _ analogous to the insulating response of 
the quantum system) . This scaling predicts that in three 
dimensions renormalized Frank constants diverge as £ 
as the crystallization transition is approached from the 
columnar phase. The singular dependence of elastic pa- 
rameters for columnar systems undergoing a continuous 
crystallization transition is summarized in Figure O 

In summary, we carried out a first-order e expansion 
below 4 dimensions to study the critical properties of the 
hexagonal-columnar liquid crystal to hexagonal crystal 



phase transition. Depending on the in-plane elastic re- 
sponse of the system, this transition proceeds either as 
first order (driven by fluctuations) or second order de- 
scribed by the critical exponents of the XY model in the 
same dimension. Finally, we argued that this rare exam- 
ple of a continuous crystallization transition, leads to a 
generic, singular evolution of bulk elastic properties. 
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